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ABSTRACT 


The  equation  of  conservation  of  mass  and  momentum  for  two-phase  flows 
are  derived  by  applying  a  generic  averaging  process.  The  properties  of  this 
averaging  process  are  discussed.  Constitutive  equations  are  proposed  for  the 

i 

inter facial  force.  In  addition ,  simple  assumptions  are  made  for  the  Reynolds 
stresses  and  the  fluid  viscosity.  These  assumptions  are  examined  for  transition 
layers . 
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CONTINUUM  MODELING  OF  TWO-PHASE  FLOWS 


Donald  A*  Drew* 


Introduction 

The  flow  of  two  materials,  one  dispersed  throughout  the  other,  has 
received  much  attention  in  recent  times*  Unfortunately,  at  this  time,  there 
seems  to  be  no  set  of  equations  which  is  regarded  as  fundamental,  from  which 
other  models  can  be  derived  as  approximations*  (Consider  the  analogy  with 
fluid  mechanics,  where  the  incompressible,  inviscid  equations  are  thought  to 
be  valid  approximations  outside  of  shear  layers  and  boundary  layers,  when 
thermal  and  sonic  effects  are  unimportant.) 

Many  researchers  derive  equations  of  motion  by  applying  an  averaging 
process  to  the  microscopic  -equations  of  motion.  The  choice  of  averaging 
process  is  dictated  by  the  taste  of  the  researcher  as  well  as  the  particular 
problem  studied.  In  this  paper,  we  give  a  derivation  of  the  averaged 
equations  by  applying  a  generic  average.  The  relation  of  the  generic  average 
to  time-  and  space-averaging  is  discussed. 

Once  believable  equations  of  motions  have  been  formulated,  it  is  natural 
to  study  their  predictions  in  relatively  simple  flow  situations.  Often  the 
constitutive  assumptions  used  in  the  model  are  derived  and/or  tested  on 
uniform  flow  situations.  We  give  a  discussion  of  transition  layers  in  two- 
phase  flows.  A  transition  layer  is  a  thin  region  where  the  concentration  of 
one  material  changes  rapidly  in  space.  An  example  is  the  "interface"  between 
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a  carbonated  liquid  and  its  foamy  "head".  An  adequate  description  of  these 
transitions  provides  a  harsh  test  of  the  constitutive  assumptions  used  in  the 
model. 

EQUATIONS  OF  MOTION 

Each  material  is  assumed  to  be  a  continuum,  governed  by  the  partial 
differential  equations  of  continuum  mechanics.  The  materials  are  separated  by 
an  interface,  which  is  a  surface.  At  the  interface,  jump  conditions  express 
the  conditions  of  conservation  of  mass  and  momemtum. 

The  equations  of  motion  for  each  phase  are  (Truesdell  and  Toupin  1960) 

( 1 )  conservation  of  mass 

■fj+V-pr-O  (1) 

(2)  conservation  of  linear  momentum 

+  V«pw  -  V*T  ♦  pf  (2) 

valid  in  the  interior  of  each  phase.  Here  p  denotes  the  density,  ▼  the 
velocity,  T  the  stress  tensor,  and  t  the  body  force  density.  Conservation 
of  angular  momentum  becomes  T  -  T*,  where  t  denotes  the  transposed.  At 
the  interface,  the  jump  conditions  are 

( 1 )  jump  condition  for  mass 

Ip(v  -  ▼i)  •*»!  -  0  (3) 

(2)  jump  condition  for  momentum 

Ipv(w  -  ▼.)••  -  T»n]  -  0  .  (4) 

Here  (  ]  denotes  the  jump  across  the  interface,  v^  is  the  velocity  of  the 

interface,  and  m  is  the  unit  normal  (Aris  1962).  we  shall  assume  that  m 

points  out  of  phase  k,  and  that  the  jump  between  f  in  phase  k  and  f  in 

Jt  k 

phase  t  defined  by  [f J  “  f  -  f  ,  where  a  superscript 

k  denotes  the  limiting  value  from  the  phase  k  side.  As  a  sign  convention 


for  the  curvature #  we  assume  that  <  is  positive  (concave)  toward  -n.  The 
mass  of  the  interface  and  surface  stresses  have  been  neglected*  We  do  not 
discuss  any  thermodynamic  relations  in  this  paper* 

Constitutive  equations  must  be  supplied  to  describe  the  behavior  of  each 
material  involved*  For  example#  if  one  material  is  an  incompressible  liquid, 
then  specifying  the  value  of  p#  and  assuming  T  *  -pi  +  U(  Fv  +  (Vv)t) 
determines  the  nature  of  the  behavior  of  the  fluid  in  that  phase*  similar 
considerations  are  possible  for  solid  particles  or  a  gas*  The  resulting 
differential  equations#  along  with  the  jump  conditions#  provide  a  fundamental 
description  of  the  detailed  or  exact  flow* 

Usually#  however#  the  details  of  the  flow  are  not  required*  For  most 
purposes  of  equipment  or  process  design#  averaged,  or  macroscopic  flow 
information  is  sufficient*  Fluctuations,  or  details  in  the  flow  must  be 
resolved  only  to  the  extent  that  they  affect  the  mean  flow  (like  the  Reynolds 
stresses  affect  the  mean  flow  in  a  turbulent  flow)* 

Averaging 

In  order  to  obtain  equations  which  do  not  contain  the  details  of  the 
flow#  it  has  become  customary  to  apply  some  sort  of  averaging  process*  We 
present  a  generic  averaging  method#  and  its  results* 

Let  <  >  denote  an  averaging  process  so  that  if  f(x#t)  is  an  exact 
microscopic  field#  then  <f>(x#t)  is  the  corresponding  averaged  field* 

An  averaging  process  assigns  average  values  to  certain  variables.  The 
ensemble#  or  set  of  possible  outcomes#  can  be  taken  to  be  the  possible  flows 
in  some  apparatus  where  the  initial  and  boundary  conditions  which  are 
prescribed  are  equivalent  in  some  sense*  For  example#  for  spherical 
particles#  it  may  be  necessary  to  give  the  statistical  distribution  of  the 
positions  and  velocities  of  the  centers  of  the  particles  at  time  t  *  0  such 
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that  the  average  number  danaity  and  average  particle  velocity  is  the  sane  for 
all  equivalent  flows.  We  shall  assune  that  there  is  some  ensemble  Q,  with 
sone  appropriate  weighting  |i(w)dw  so  that  the  average  of  f  is  given  by 

<f>(n,t)  -  /Q  e(K,t,u)tU  w)du  . 

Two  cases  can  be  discussed.  If  the  flow  is  nearly  steady,  so  that  a  time 
translation  T  makes  no  essential  difference  in  the  ensemble,  it  may  be 
enough  to  consider  the  subset  of  the  entire  ensemble  which  consists  of 
translations  in  time  of  amount  T.  We  assign  a  weight  m(  t)  to  the 
likelihood  of  the  flow  whose  outcome  at  time  t  is  f(*,t-  t),  where 
f(x,t)  is  the  outcome  at  time  t  in  some  flow.  The  average  of  f  is  then 
taken  to  be 

<f>t(*,t)  -  f<m,t-T)M(  T)dr  . 

This  is  classical  tiM  averaging;  it  is  often  used  with 

f i  if  0  «  *  <  * 

y(  t)  «  \ 

V  0  otherwise  , 

although  other  averages  are  peesible. 

If  there  are  no  boundaries  in  the  flow  (that  is,  boundary  effects  are 

unimportant ) ,  then  small  spatial  translations  should  make  no  difference  in  the 

ensemble.  In  analogy  to  the  above,  the  average 

<f>  (*#t)  ■  /  3  f(*+m,t)i*(m)ds 
R 

can  be  defined.  This  is  the  classical  space  average;  it  is  often  used  with 

fi  “  •«* 

m(s)  -  \ 

\  0  otherwise 
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where  v  la  son*  volume  (for  example,  a  sphere).  Again,  other  averages  are 
possible. 

Thus,  in  scow  sense,  ensemble  averaging  contains  space  and  time  averages 
as  special  cases.  The  averaging  process  is  assumed  to  satisfy 


<f+g>  “  <f>  +  <g> 

(5) 

«f>g»  -  <f><g> 

(6) 

<c>  ■  c 

(7) 

it  a 
"  "at 

(8) 

it  i 

<f>  • 

(9) 

The  first  three  of  these  relations  are  called  Reynolds  rules,  the  fourth 
is  called  Leibnitz'  rule,  and  the  fifth  is  called  Gauss'  rule. 

In  order  to  apply  the  average  to  the  equations  of  motion  for  each  phase, 
we  introduce  the  phase  function 
Xj^Xft)  which  is  defined  to  be 

f  1  if  x  is  in  phase  k  at  time  t 
X^(x»t )  -  <  (10) 

\  0  otherwise  . 


We  shall  deal  with  X^  as  a  generalized  function,  in  particular  in  regard  to 
differentiating  it.  Recall  that  a  derivative  of  a  generalized  function  can  be 
defined  in  terms  of  a  set  of  "test  functions"  +,  which  are  "sufficiently 


3Xk 

smooth"  and  have  compact  support.  Then  and  -53-  are  defined  by 


/ 


K 

ax. 


R3XR 


■j~  (x,t)  +(x,t)dxdt  - 


-  -  /  ,  X  (x,t) 
R3W  * 


1* 

at 


(x,t)dxdt  , 


(11) 
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3X. 

/  ,  (*»t) $(x,t)dxdt  - 

RXR 


-  /  3  X^Xft)  -|£-  (x,t)dxdt  . 


It  can  be  shown  that 


■§r+  v*k" 0 

in  the  sense  of  generalized  functions. 

If  f  is  smooth  except  at  S,  then  f  VX^  defined  via 


/  3  fVX^^dxdt  -  £m  fs  i^fN  dSdt  , 


where  is  the  unit  normal  exterior  to  phase  k,  and  f*  denotes  the 

limiting  value  of  f  on  the  phase-k  side  of  s. 

It  is  also  clear  that  Vx^  is  zero,  except  at  the  interface.  Equation 
(14)  describes  the  behavior  of  VX^  at  the  interface.  Note  that  it  behaves 
as  a  "delta-function",  picking  out  the  interface  S,  and  has  the  direction  of 
the  normal  interior  to  phase  k. 

Averaged  Equations 

In  order  to  derive  averaged  equations  for  the  motion  of  each  phase,  we 
multiply  the  equation  of  conservation  of  mass 
valid  in  phase  k  (1)  by  X^  and  average.  Noting  that 

3k 

'Seif  “*kv- p  dT“-kv+  (15) 


we  have 


\  V»pv  -  -  pr*VXk  , 


*§£  <\P>  +  “  <[p(^Ti)]k*VXk> 


Similar  considerations  for  the  momentum  equations  yield 
3 

<xjcpir>  +  V*<Xkpw>  -  +  <xk<3f> 

+  <  [  )  -  T]k‘VXk>  . 
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The  terms 


<[p(e-vi)]K*V5ck>  -  rk  <19> 

and 

<[pv(e-Ti>  -  «k*7xk>  -  (20) 

are  the  interfacial  source  terms*  As  noted ,  picks  out  the  interface, 

and  causes  discontinuous  quantities  multiplying  it  to  be  evalued  on  the  phase- 
k  side  of  the  interface. 

The  jump  conditions  come  from  equations  (3)  and  (4).  We  have 
2  v  2 

l  <[p(w.  )]K»7)t  >  -  i  r  •  0  (21) 

k-1  Jp-1 


£  <(pv(v-v^)  -  *]  •V3Cjt>  ■  0  .  (22) 

k*1 

Applying  a  more  specific  averaging  process  (time  averaging,  for  example) 
requires  a  different  set  of  manipulations  regarding  the  interfacial  source 
terms  (Anderson  &  Jackson  1967,  Drew  1971,  Ishii  1975,  Delhaye  and  Achard 
1976).  Almost  all  of  the  derivations  for  specific  averaging  processes  seem  to 
be  more  complicated  than  the  above;  however,  the  trade-off  for  the  simple 
derivation  is  that  all  manipulations  now  involve  generalised  functions. 

The  volumetric  concentration  (or  volume  fraction,  or  relative  residence 


time)  of  phase  k  is  defined  by 


We  note  that 


\  ■  ‘V 


»t  »t 


•  ‘V  • 
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There  arc  two  typos  of  averaged  variables  which  are  useful  in  two-phase 
mechanics #  namely  the  phasic#  or  X^-weighted  average#  and  the  mass-weighted 
average*  Which  is  appropriate  is  suggested  by  the  appearance  of  the  quantity 
in  the  equation  of  motion*  The  phasic  average  of  the  variable  +  is  defined 
by 

\  -  <«> 

and  the  aaw  weighted  average  of  the  variable  ,  is  defined  by 

\  *  •  <27> 
It  is  convenient  to  write  the  stresses  in  terms  of  pressures  plus 

extra  stresses e  Thus, 

\  -  -Pkl  +  ^  .  (28) 

It  is  expected  that  readers  familiar  with  fluid  dynamical  concepts  are 
familiar  with  the  concept  of  pressure  in  fluids;  in  this  case,  p^  can  be 
thought  of  as  the  average  of  the  microscopic  pressure  *  If  one  of  the  phases 
consists  of  solid  particles#  the  concept  is  less  familiar*  In  this  case,  the 
microscopic  stress  (involving  small  elastic  deformations,  for  exas^le)  is 
thought  of  being  made  up  of  a  spherical  part  (actinj  equally  in  all 
directions)  plus  an  extra  stress*  The  spherical  part,  when  averaged,  yields 
the  pressure  p^  in  equation  (28)* 

It  has  further  become  customary  to  separte  various  parts  of  the 
inter facial  momentum  transfer  term*  This  is  done  by  defining  the  inter facial 
velocity  of  the  phase  by 

W,i  ■  .  <»»> 

and  tha  lntarfacial  pressure  on  the  k*®*  phase  by 

pk(i,,\|2  ■  V’\  •  <30) 

■quatlon  (30)  is  the  dot  product  of  Vc^  of  the  "standard"  definition  (I3hii 
1975)  of  the  interfacial  pressure.  The  standard  definition  uses  three 
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aquations  to  define  one  scalar  quantity*  and  cannot  be  a  generally  valid 
definition.  Here  the  remaining  part  of  the  contribution  of  the  pressure  at 
the  interface  is  lumped  with  the  viscous  stress  contribution  at  the  interface* 
and  is  treated  through  the  use  of  a  constitutive  equation.  Thus,  we  write 

^  “  Vk,i  "  P*,iV<V  +  "fc  «  (31) 

d  )c  u 

where  «  <(p-pfc  *VX^>  is  referred  to  as  the  interfacial  force 

density*  although  it  does  not  contain  the  effect  of  the  average  force  on  the 
interface  due  to  the  average  interfacial  pressure.  The  term  -p^  ^Vo^,  which 
does  contain  the  force  due  to  the  average  interfacial  pressure,  is  sometimes 
referred  to  as  the  bouyant  force.  The  reason  for  this  terminology  is*  of 
course*  that  the  buoyant  force  on  an  object  is  due  to  the  distribution  of  the 
pressure  of  the  surrounding  fluid  on  its  boundary. 

With  equations  (23)  and  (26)  -  (31),  the  equations  of  motion  (17)  and 


(18)  become 


*  ''w  -  rk  <“> 


‘Wk 


* 'm  ■  ♦  »*v  \  *  v 


♦  rk\.i  *  "’k.i  -  »k,7\ 


*4- 


The  jump  conditions  (21)  and  (22)  are 


l  r.  -  o 


i,  irk\,i  *  “k.iH  *  "k1 


o 


(35) 


Adequate  models  for  compressibility  and  phase  change  require 
consideration  of  thermodynamic  processes*  These  are  beyond  the  scope  of  this 
paper;  therefore  we  shall  restrict  our  attention  to  incompressible  materials 
where  no  phase  change  occurs*  Thus  we  assume  that 

*  constant  (36) 

and 


rk  *  0  .  (37) 

In  order  to  simplify  the  notation,  we  shall  drop  all  symbols  denoting 
averaging* 

CONSTITUTIVE  EQUATIONS 


In  order  to  have  a  useable  model,  relations  must  be  given  which  specify 
the  stresses  (x  +  c^),  the  interfacial  force  density  M^,  and  the  pressure 
differences  p^  -  consistent  with  the  equations  of  motion  and  the  jump 

conditions* 

The  fundamental  process  consists  of  proposing  forms  for  the  necessary 
terms  within  the  framework  of  the  principles  of  constitutive  equations, 
finding  solutions  of  th  *  resulting  equations,  and  verifying  against 
experiments.  The  ideal  end  result  of  the  process  is  a  set  of  equations  which 
could  be  used  to  predict  the  behavior  of  the  two-phase  flow,  for  example  with 
a  computer  code*  With  the  equations  should  come  a  set  of  conditions  for  the 
validity  of  the  values  of  the  constants  and  other  functions  used  in  the 
constitutive  equations* 

The  stresses  the  inter  facial  force  density 

d 

and  the  pressure  differences  Pk  "  ^k  i  are  *8sumed  to  be  functions  of 


o^,  dc^/dt,  Vc^,  ▼k,  7*k,  3*k/3t,  ...  , 


where 


•  •  • 


represents  the  material 


properties,  such  as  the  viscosities  and  densities  of  the  two  materials,  and 
other  geometric  parameters  such  as  the  average  particle  sire,  or  the 
interfacial  area  density* 

For  concreteness,  we  shall  refer  to  phase  one  as  the  particulate,  or 
dispersed  phase  and  include  in  that  description  solid  particles,  droplets,  or 
bubbles*  Phase  two  is  then  the  continuous,  or  carrier  phase,  and  can  be 
liquid  or  gas*  We  shall  denote 

a  «  ot  (38) 

so  that 

1-a  -  •  (39) 

It  is  evident  that  both  a  and  1-a  need  not  be  included  as  independent 
variables  in  forming  constitutive  equations* 

Drew  and  Ijahey  (1979)  consider  the  general  process  of  constructing 
constitutive  equations*  The  simplest  reasonable  set  of  assumptions  leads  to 
models  for  the  motions  of  the  two  materials  which  may  be  ill-posed*  Drew 
(1982)  gives  a  review  of  the  state  of  affairs* 

Essentially,  the  problem  and  the  reason  for  its  importance  can  be 
summarized  in  the  following  manner*  The  simplest  model  assumes  that  the 
interaction  forces  are  due  to  viscous  drag,  and  that  pressure  forces 
equilibrate  across  particles  instantaneously*  Thus, 

Q 

M1  *  "  8  *P1  'T  'VT2|(  W  (40) 

and 


*1  “  *2  “  *1,i  “  *>2,i  *  P  •  (41) 

If  the  stresses  are  ignored,  the  model  has  complex 

characteristics,  and  hence  is  ill-posed*  see  Ramshaw  and  Trapp  (1978).  The 
ill-posed  nature  of  the  model  leads  (theoretically)  to  solutions  which  grow 
rapidly  on  a  small  length  scale*  This  leads  some  to  conjecture  that  the  model 


-11- 


is  trying  to  resolve  events  on  the  microscale.  The  application  of  an  average , 
however,  is  supposed  to  lose  details  of  the  flow;  indeed,  both  (40)  and  (41) 
specifically  ignore  flow  details  on  the  scale  of  the  particle  size  and 
smaller.  The  obvious  conclusion  is  that  equations  (32  -  37),  with  assumptions 
(40),  (41)  and  neglect  of  viscous  forces  lack  some  mechanism  which  is 
important  on  the  particle  scale.  It  is  noteworthy  that  no  difficulties  arise 
in  numerical  simulations  with  the  ill-posed  system.  The  reason  for  this  is 
that  the  instability  inherent  in  the  ill-posed  system  appears  on  a  scale 
comparable  to  the  particle  radius  to  a  power  which  depends  on  the  exact  form 
of  the  drag  law  used  for  Cp.  Thus,  if  the  mesh  is  not  refined  to  the 
particle  size,  no  instabilities  will  be  seen.  This  suggests  that  computation 
with  the  ill-posed  system  with  a  reasonable  mesh  will  most  likely  give  good 
results. 

The  above  argument  indicates  some  difficulties  for  stability 
calculations.  When  is  an  instability  not  an  instability?  Presumably  it  is 
unobserved  if  its  wavelength  is  too  small.  It  would  be  satisfying  to  find  the 
missing  effect  in  the  model  and  show  how  the  model  reduces  to  the  isobaric, 
inviscid  model  mentioned  above.  There  is  a  long  list  of  candidates  for  the 
forgotten  mechanism  (Drew  1982)  but  no  clear  winner  has  emerged.  The  sensible 
approach  seems  to  be  to  examine  various  models  on  a  mesoscale,  that  is,  on  a 
scale  which  is  small  compared  to  the  usual  experimental  verification  flows 
(such  as  laminar  settling  in  a  still  fluid),  but  large  compared  to  the 
particle  scale.  These  mesoscale  calculations  may  indicate  what  sorts  of  terms 
are  needed  to  give  a  valid  description  on  a  smaller  scale. 

There  is  another  reason  for  seeking  the  more  complete  model.  If  the  ill- 
posed  model  is  the  limit  of  some  more  complete  model  with  some  effects 
neglected,  if  the  neglect  can  be  done  by  a  set  of  formal  manipulations  on  the 
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aore  complete  model,  the  result  may  Indicate  whether  some  well-posed  model 
might  do  as  well  as  the  above  Ill-posed  model.  Moreover,  it  is  always  of 
interest  to  find  reduced  models  which  contain  the  same  features  as  the 
original  model. 

He  shall  examine  the  effects  of  the  viscous  and  Reynolds  stresses  on  a 
mesoscale  motion,  namely  transition  layers  in  vertical  flow.  If  the  flow  is 
vertical,  and  we  denote  the  vertical  velocities  by 

▼l  ■  u(s,t)k  (42) 

▼2  ■  v(*,t)k  ,  (43) 


we  have 


8a  8ou 
8t  8* 

au  +  (1-a )v 


0 

j(t) 


(44) 

(45) 


/in  ^  8u\  8d  .  . 

“V*  +  u  “  "a  fe  +  0b(VWu)  " 


(46) 


"  +  fc  °°1 


(1-a)p2(|^  +  v  -|-»)  -  -(1-a)  -j£  +  ob(u-v)  -  (1-a)p2 g  + 


♦  t  ^  -5  • 


(47) 


The  function  j(t)  is  the  volumetric  flux,  and  will  be  taken  to  be 
constant.  The  stress  models  are 

»  0  (48) 

TJ  ■  •  -5  ««> 

to  represent  the  viscosity  of  the  fluid.  The  particles  are  assumed  to  be 
inviscid.  In  addition,  the  Reynolds  stress  terms  will  be  taken  to  be 

constants. 
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For  problems  of  sedimentation  and  transitions  in  fluidized  beds,  we  shall 


assume  that  j(t)  -  constant  for  t  >  0. 

Kinematic  Waves 

If  we  ignore  inertia  in  equations  (46)  and  (47)  and  set  =  0  and 

M  *  0,  we  have  a  model  which  reduces  to  a  one-dimensional  scalar 
conservation  law.  If  the  pressure  is  eliminated  from  the  resulting  momentum 
equations,  we  have 

( 1  -  a)  ( p  -  p  )  g 

v— u  -  - ^ — —  .  (50) 

Using  (45)  gives 

u  -  j  -  (1-a)(v-u)  • 


Hence 


au 


aj  - 


o(1-a)  (P1-P2)g 


f  ( a) 


(51) 


Equation  (44)  with  au  given  by  (51)  is  a  scalar  conservation  law  for 


a(z,t). 

Solutions  can  be  found  by  the  method  of  characteristics;  on  the 
characteristics. 


da 

dt " 0 


(52) 


dz 


—  ■  f'(a)  *  constant  . 
at 


(53) 


As  long  as  no  discontinuities  occur,  a  is  constant  on  characteristics. 


If  discontinuities  occur,  they  must  propagate  at  a  speed  given  by 


la)  ' 


(54) 


where  [ 4)  denotes  the  jump  in  the  quantity  +  across  the  discontinuity.  A 
discontinuity,  or  shock  will  persist  if  characteristics  tend  to  come  into  the 
shock.  If  character-  istics  diverge  from  the  shock,  then  it  will  smooth 
out.  in  Figure  1,  a  shock  from  (a  ,f(a  ))  to  (a  ,f(a  ))  is  not  stable  if 


the  region  where  occurs  is  above  the  shock*  This  is  because  the 

characteristics  at  are  moving  faster  than  the  shock*  Note  that  contact 

discontinuities  are  possible  with  this  model* 


Figure  I 

For  fluidization,  j  >  0.  Zn  order  to  describe  a  transition  in 
fluidization,  we  assume  that  the  bed  is  fluidized  at  t  *  0-  with  a 
concentratl  a.  o+  which  corresponds  to  some  value  of  the  volumetric  flow 
rate  j+.  At  time  t  =  0,  the  volumetric  flow  rate  is  instantaneously 
changed  to  j.  The  flow-concentration  diagram  is  given  in  Figure  1.  The 
simplest  situation  is  depicted,  consisting  of  an  upward  traveling  bottom  shock 
(transition)  and  an  upward  traveling  top  shock. 

Figure  2  shows  the  consentration  at  some  time  t1 ,  and  Figure  3  shows 
the  solution  in  various  regions  of  the  t  -  z  plane. 
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Figure  3 


Dlffuaional  Regularization 

One  way  to  understand  shocks  is  to  include  in  the  model  some  means  of 
smearing  them  out.  This  can  be  done  in  the  present  problem  by  Including  the 
diffusion  terms.  We  continue  to  ignore  inertia  and  viscosity.  Repeating  the 
procedure  outlined  above  yields 

a»  -  f(o)  -  D  -g  ,  (55) 


where 


(56) 


(l-a)l(1-o)o1  +  aa2J 


3a  Hi  a)  3_ 

at  a*  "  at  ^  a* 


The  equation  for  a  is  a  nonlinear  diffusion  equation 

•£  (■>  £)  •  <5” 
In  order  to  examine  the  transition,  we  let 

o  -  aiz  -  st)  (58) 

where  s  is  the  speed  of  propagation  of  the  traveling  wave  representing  the 

transition.  We  obtain 

Sr  tf(o>  -  .<■)  -jh  .  (59) 


Integrating  from  -•  to  £  gives 

f  la)  -  [f  (a_)  +  s(  a  -  aj]  »  D  — jp  .  (60) 

The  quantity  in  the  bracket  is  the  equation  for  the  chord  (in  Figure  1  for 
example) ,  and  f(a)  is  the  flow-concentration  curve  there*  If  the  curve  lies 
above  the  chord ,  then  do/d£  is  positive*  If  the  curve  lies  below  the  chord, 
then  dq/d£  is  negative*  The  transitions  go  from  a_  to  a+  where  is 

the  value  of  a  where  the  left  hand  side  of  (60)  is  zero,  that  is,  where  the 
curve  and  chord  intersect*  If  D  is  small,  the  transition  region  is  thin* 

The  diffusional  regularization  corresponds  to  the  results  obtained  from  shock 
stability  considerations* 

Inclusion  of  the  inertia  of  both  phases  complicates  the  situation 
immensely*  It  can  be  shown  that  the  thing  which  corresponds  to  the  shock  in 
the  kinematic  wave  model  is  not  a  shock  in  the  model  which  includes  inertia* 

In  order  to  study  the  transition,  assume  that  a,  u,  v  and  p  are 
functions  of  £  *  z  -  st,  where  s  is  the  speed  of  the  transition  wave,  and 
o  ♦  o_,  u  ♦  u_,  v  ♦  v_  as  £  ♦  •,  o  ♦  a+,  u  ♦  u^,  v  ♦  v^  as  £  ♦  •* 

V  u^  and  y±  are  related  to  s  by  (45)  and  (54).  Substituting  in  the 
equations  (44  -  47)  and  eliminating  p,  u,  and  v  gives 
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m(i-o  )(v  -s) - -  or  + 


(1-a) 


(61) 


+  B(  a)  o'  +  g(  a)  -  0 


where 


8(a) 


a  °-  p1(u-’s)2 

0+0  -  -  » 

1  2  1-a  2 


a  t e  v2,  x2 

p.(1-a  )  (v  -s) 


(1-a) 


3  '2' 


(62a) 


and 

g(a)  - ~ — -  [f  (  a)  -  a  u  -s(  a- a  )]  .  (62b) 

(1-a) 

Note  that  if  u  33  0,  and  we  wish  to  retain  the  picture  given  by  the 
characteristics  of  the  scalar  conservation  law,  we  can  do  so  only  if 

A 

8(a)  >  0  (63) 

for  a  between  a_  and  a+.  In  a  sense,  equation  (63)  suggests  that  inertia 
will  be  unimportant  for  transitions  If  diffusion  Is  sufficiently  large; 
sufficiently  large  means  so  that  equation  (63)  is  satisfied* 

Note  from  equation  (61)  with  p  »  0  and  *  0,  that  inertia 

without  diffusion  gives  results  opposite  to  the  results  from  the  scalar 
conservation  law* 

The  Effect  of  Viscosity 

dw 

In  equation  (61),  let  a*  -  w(a)*  Then  *  w  — ,  and  equation  (61) 

aft 

becomes 


p(1-a  ) (v  -s)  — 2 — ~  +  M(1-a  ) ( v  -a)  , 

‘  ‘  (1-a) 2  d°  ‘  ‘  (1-a)3 


a  2 
w 


(64) 


+  B(  a)w  +  g(  a)  ■  0  . 


Let  us  now  define  G(  a)  by 


and 


«. 


G(a_)  »  0 

if 

g(a)  >  0 

for 

a_  <  a  <  a^ 

G(a+)  »  0 

if 

g(  a)  <  0 

for 

a_  <  a  <  a^ 

(We  also  assume  that  (1-a)(v^s)  >  0.)  Equation  (64)  becomes 


a  i  2 

dS  [2  ^  +  G(a)J  “  *6(<I,W  ‘  (67) 

1  w  2 

The  curves  H(  a,w)  -  -(- — )  +  G(  a)  are  closed  curves  centered  at  a  if 

2vi-cr  • 

if  q(a)  >  0,  and  a+  if  g( a)  <  0.  See  Figure  4.  On  a  trajectory  leaving 
a+  (if  g(a)  >  0)  or  a_  (if  g( a)  <  0),  the  function  a(  £)  satisfies 

H(a,w)  -  -0(a)w2  (68) 

which  is  negative  if  8(a)  >  0.  Thus  H  decreases  on  a  trajectory,  giving 
the  transition  from  to  a_  (if  g(  a)  >  0)  or  from  a_  to  (if 

g(a)  <  0 )  *  Since 


g(a)  - - - -  (f  (  a)-a_u_  -  s(  »-a_) )  , 

a(  1-a#  ,  y(1-a_)  (v_-s) 

we  see  that  if  the  curve  f  is  above  the  chord,  s(a-a_)  +  a_u_,  then 
g(a)  <  0,  and  this  type  of  transition  can  occur  from  a_  to  a+*  On  the 
other  hand,  if  the  curve  f(a)  is  below  the  chord  s(a»a_)  +  a_u_,  then 
g(a)  >  0,  and  the  transition  can  occur  from  a+  to  a_.  This  again  agrees 
with  the  picture  given  by  the  characteristics  for  the  one-dimensional 
conservation  law* 

Conclusion 

There  are  two  distinct  features  of  two-phase  flow  modeling  addressed  in 
this  paper,  namely  the  averaging  process,  and  transition  layers* 

The  generic  averaging  process  includes  certain  classical  averages  as 
special  cases  and  gives  the  same  results  as  these  classical  averages  when 
applied  to  the  equations  of  motion* 
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Figure  4 


Transition  layers  are  often  observed  in  two-phase  flows,  and  various 
models  for  them  have  been  used.  The  simplest  model,  that  of  the  kinematic 
wave,  describes  the  situation  when  inertia,  viscosity  and  diffusivity  are 
negligible,  but  diffusivity  dominates  over  inertia. 


REFERENCES 


Anderson,  T.  B.  and  Jackson,  R.  1967*  A  fluid  mechanical  description  of 
fluidized  beds,  Ind.  Eng*  Chem.  Fund*  6,  527-539, 

Aris,  R.  1962.  Vectors,  Tensors  and  the  Basic  Equations  of  Fluid  Mechanics, 
Prentice  Hall,  Englewood  Cliffs,  NJ. 

Delhaye,  J.  M.  and  Achard,  J.  !*•  1976.  On  the  averaging  operators  introduced 
in  two-phase  flow  modeling,  CSNI  Specialists  Meeting  on  Transient  Two- 
Phase  Flow,  Toronto,  Canada. 

Drew,  D.  A*  1971*  Average  field  equations  for  two-phase  media.  Studies  in 
Appl*  Math*  L: 133-166. 

Drew,  D.  A*  and  Lahey,  R*  T*  1979.  Application  of  general  constitutive 
principles  to  the  derivation  of  multi-dimensional  two-phase  flow 
equations,  Int*  J*  Multi-phase  Flow,  5:243-264. 

Drew,  D.  A*  1982.  Mathematical  Modeling  of  Two-Phase  Flow,  to  appear  in 
Annual  Review  of  Fluid  Mechanics. 

Ishii,  M.  1975*  Thermo-fluid  Dynamic  Theory  of  Two-Phase  Flow,  Eyrolles, 
France • 

Ramshaw,  J.  D*  and  Trapp,  J.  A*  1978*  Characteristics,  stability  and  short- 
wavelength  phenomena  in  two-phase  flow  equation  systems.  Nuclear  Sci. 
and  Eng*  66:93-102. 


DAD/ j vs 


21- 


SSCUMYV  CkMHnCtTNN  OF  THIS  PASS  &Ma  Mm 0 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


.  RSCIRISNT'S  CATALOG  NUMBER 


•2438 


«.  riTLK  MMt) 

Continuum  Modeling  of  Two-Phase  Flows 


I  fcM  TTmrrr 


Donald  A.  Drew 

i.  PERFORMING  ORGANIZATION  NAME  ANO  AOORES8 

Mathematics  Research  Center,  University  of 

610  Walnut  Street 

Wisconsin 

Madison.  Wisconsin  53706 

II.  CONTROLLING  OFFICE  HANK  AND  ADDNKSS 

U.  S.  Army  Research  Office 

P.O.  Box  12211 

Research  Trlanole  Park.  North  Carolina 

27709 

i.  performing  org.  report  number 


8.  contract  or  grant  numbers 
DAAG29-80-C-0041 


12.  report  date 
October  1982 


«2.  NUMBER  OP  PAGES 
21 


NAME  *  ADDREtSfff  different  fra  Controlling  OtHco)  IS.  SECURITY  CLASS,  (of  (fife  report; 

UNCLASSIFIED 


ASSIFICATION/OOWNGRADING 

OULE 


«.  DISTRIBUTION  STATEMENT  (of  this  Repeat) 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  die  obotroci  entered  In  Block  SO.  It  different  ft cm  Report) 


IS.  KEY  WORDS  (Contlnuo  on  reeorae  cldc  It  noceaaery  and  Identity  by  bloc *  number) 

two-phase  flows,  transition  layers,  shocks,  modeling 


20.  ABSTRACT  (Continue  on  rarer**  aide  It  noceaaery  and  identity  by  block  mmbat) 

The  equation  of  conservation  of  mass  and  momentum  for  two-phase  flows 
are  derived  by  applying  a  generic  averaging  process.  The  properties  of  this 
averaging  process  are  discussed.  Constitutive  equations  are  proposed  for  the 
interfacial  force,  in  addition,  simple  assumptions  are  made  for  the  Reynolds 
stresses  and  the  fluid  viscosity.  These  assumptions  are  examined  for  transitior 
layers. 


I  JAM  71 


COITION  OP  I  NOV  M  It  OStOLtTC 


UNCLASSIFIED 

MCUNITY  CLASSIFICATION  OP  THIS  MAO 


